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Nonequilibrium dynamics at interfaces is generally driven by a chemical potential. Here we 
demonstrate a generic technique to derive the basic equations of motion, boundary conditions and 
the chemical potential in a consistent way from fundamental variational principles. As a particular 
example, we consider a solid surface with elastodynamic effects, together with surface energy and 
tension. We apply the generic results to perform a linear stability analysis of a planar front subjected 
to uniaxial stress (Asaro-Tiller-Grinfeld instability Q), here also with surface tension. 
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I. INTRODUCTION 

The quest for the correct equations of motion describ- 
ing a system is an essential task of each problem in all 
branches of physics. The procedure is well-known for 
mechanical systems, where for example the appropriate 
Lagrangian has to be found, from which the (Newtonian) 
equations of motion can be derived using variational prin- 
ciples. This method can in particular be applied to con- 
tinuum mechanics, generating the usual bulk equations 
and boundary conditions. 

However, for more complicated dissipative systems, 
these straightforward and general methods are usually 
not applicable. Bulk diffusion, friction, and also irre- 
versible phase transitions are important examples. In the 
framework of linear nonequilibrium thermodynamics, the 
"process velocity" is proportional to the deviation of an 
energy functional from its equilibrium value. 

If we consider for example diffusion processes, then the 
flux of material is typically j = — DV/z with the diffu- 
sion coefficient D and a chemical potential /i. Similarly, 
for phase transformations at an interface, the resulting 
interface normal velocity is assumed to be proportional 
to the difference in chemical potentials in both phases, 
Vn = wA/i (k is a kinetic coefficient); obviously, in both 
cases the chemical potential /i is the relevant driving force 
for this type of dynamics; however, it is not evident how 
it can be expressed through the state variables of the 
system. Therefore it is desirable to provide a method 
to derive the correct expression for the chemical poten- 
tial using flrst principle techniques. For relatively simple 
processes, it is easy to guess it correctly, but as soon 
as more complicated effects are included, the results are 
sometimes unexpected. A particular example are elasto- 
dynamic effects, and the correct expression for the chem- 
ical potential will be derived in this paper. The influence 
of surface tension, also causing complicated contributions 
to the chemical potential, was already considered for pro- 
cesses at elastically stressed surfaces in 0; however, the 
author made some systematic mistakes. This illustrates 
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that the lack of a consistent derivation and the use of bare 
intuition can easily lead to wrong predictions. Though 
the aim of the current paper is mainly to demonstrate 
methodological aspects how to derive consistently these 
expressions from variational principles, we also want to 
point out that the results form a basis to the understand- 
ing of important topics like solidification of superfluid 
helium and fracture mechanics Q- 

The Asaro-Tiller-Grinfeld (ATG) instability is a par- 
ticular example of these processes [ij . Here one considers 
a solid piece of material which is either in contact with 
its melt such that phase transitions and thus motion of 
the interface are possible, or which suffers from surface 
diffusion. A uniaxial stress at the surface of this solid 
gives rise to a morphological instability; long wave defor- 
mations lead to a reduction of the total energy. A relax- 
ation of elastic stresses in some regions makes them more 
favorable in comparison to the valleys in which the stress 
is increased; consequently, a corrugated surface is char- 
acterized by a modulation of the elastic energy density 
at the surface. This contribution favors a destabiliza- 
tion of the interface, counteracted by the surface energy 
which prefers a flat surface. After the initial develop- 
ment, the originally slow process becomes extremely fast 
and leads to an unphysical breakdown of the analytical 
and numerical predictions: Deep groves with sharp tips 
advancing with increasing velocity can form in the sam- 
ple ^5,, _6J . Here strong surface tension and also dynamical 
effects become relevant, but they are not taken into ac- 
count in conventional descriptions. The main problem is 
here that it is not obvious how these effects enter into 
the expression of the chemical potential. 

The main goal of the current paper is therefore to pro- 
vide a method to derive all describing quantities in a 
consistent way from variational principles. A particular 
illustration is the Grinfeld instability with surface ten- 
sion. The paper is organized as follows: In Section HTl we 
present the Lagrangian which is used to derive all further 
physical laws. In particular, we consider a linear elastic 
medium with dynamical effects, surface tension and sur- 
face energy, under the influence of external forces. In 
Section IIIII we derive the basic equations of motion and 
boundary conditions for the "mechanical" degrees of free- 
dom, the elastic displacements, from the Lagrangian. In 
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Section IIVI the expression for the chemical potential at 
the surface of a solid body is derived using variational 
principles. The surface tension enters into the final ex- 
pression in a complicated way, and also inertial effects 
appear unexpectedly. Therefore, in Scction^thc rate of 
change of the total energy during a reshuffling of the sur- 
face is calculated; this confirms the preceding results for 
the chemical potential. A specific application of these 
general techniques is the Grinfeld instability with sur- 
face tension fSection WJ^ . We derive the spectrum of 
the initial stage of this instability. Additionally, in Sec- 
tion the energetic change during the morphological 
transformations is calculated, correcting inconsistencies 
in 0] . Appendix ^ explains the fundamental conser- 
vation laws for momentum and angular momentum for 
elastic systems with surface tension; Appendix ^ con- 
tains the mathematical details to express the dissipation 
rate and the functional derivative for a generic functional 
which is used to derive the chemical potential and the en- 
ergy dissipation. 



II. ACTION AND LAGRANGIAN 

We consider a solid body in vacuum under the influence 
of external forces, linear elastodynamics, surface energy 
and tension. The Lagrangian is given by 

L= J iT-Ub)dV- J UsdT = Lb + Ls, (1) 



V{t) 



dV{t) 



where V{t) is the time-dependent volume of the elastic 
body and dV{t) its boundary with elements dr. The 
expressions for kinetic and potential energy densities are 
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Here p denotes the mass density, which we assume to be 
constant. Therefore we neglect compressive effects, since 
Ap = puii <C p- The elastic displacement field is de- 
noted by Ui, Gik is the stress and Uij- the strain tensor, 
a is the surface energy density and (3 the surface tension 
coefficient. The representation of the surface tension con- 
tribution is unusual, and we will further comment on it 
below 0. The stresses cr^^^ are the given external forces 
which may act on some parts of the boundary; every- 
where else they are zero. We assume the volume to be 
two-dimensional; r is the tangential and n the normal 
direction at the interface. 

The surface energy in Eq. expresses the energy of 
the original undeformed surface, whereas the surface ten- 
sion contribution gives the additional energy which arises 
from elongation of the surface due to elastic effects. No- 
tice that Urrdr is the elongation or compression of the 



length element dr. Surface energy gives only an addi- 
tional contribution to the total energy, surface tension 
also couples to the elastic fields and changes the elastic 
boundary conditions, as will be seen in the next section. 
Also, in Section IVIII it will be shown how the surface 
tension can be interpreted in a more intuitive way, by 
relating it to surface forces which are proportional to the 
local curvature of the boundary. In general, both coeffi- 
cients a and (3 are of the same order of magnitude, and 
for a liquid they can be shown to be identical ; in gen- 
eral, (3 can even be negative. For three-dimensional situ- 
ations, the surface tension contribution in Eq. Q reads 
in general Up = f3ijUij with the summation of i,j being 
restricted to the two tangential directions of the surface. 
In the case of surfaces with a symmetry axis higher than 
twofold, the surface tension tensor is diagonal, Pij = pSij 
7]. 

Stress and strain are connected by Hooke's law. For 
simplicity, we restrict it to the isotropic case. 
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with E being the elastic modulus and v the Poisson ratio. 
The inverse transformation is 



The action is given by the time integral 



Ldt = Sb + Ss 



(6) 



(7) 



with arbitrary initial and end times toi^i- 

Since the Lagrangian is translational and rotational in- 
variant in the absence of external forces, total momentum 
and angular momentum are conserved. This property is 
not obvious if nonequilibrium processes at the surface lo- 
cally change the amount of material; a further discussion 
is given in Appendix EI 



III. ELASTODYNAMIC EQUATIONS AND 
BOUNDARY CONDITIONS 



The Euler-Lagrange conditions demand that the vari- 
ation of the action with respect to the displacements has 
to vanish 0. Here we also take into account that the 
volume can change in time due to other processes. We 
obtain for the bulk contribution 
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Now we examine the first term 
tl 

SA ~ J '^^ J PUiSuidV 

to V(i) 

in more detail. We start from the expression 
d 
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where w„ is the normal velocity of the interface; it is pos- 
itive if the solid extends. We point out that the elastic 
velocity u„ at the surface and u„ are different quanti- 
ties: The latter one describes the morphological changes 
of the shape due to nonelastic effects, for example by 
surface diffusion or the growth of a solid body from its 
undercooled melt. On the other hand, the normal dis- 
placement rate u„ does not contribute to w„ , since in the 
framework of the linear theory of elasticity one does not 
distinguish between the undeformed and deformed frame 
of reference. Therefore we obtain 
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The first bulk integral does not contribute since we have 
to demand that all variations disappear at to and ti . Con- 
sequently, we obtain the variation of the bulk action 



SSb = / dt 
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{Uni + VnPU^^Ui dr 
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We define the curvature of the surface to be negative for 
a convex body. The tangential strain is 



dur dur 
dr dr 
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Here, the definitions of the partial and the total tangen- 
tial derivative are 

dur y + At) — Ut-(t) 

dr At^o At ' 
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Thus we get for the variation of the boundary action 
tl 

5Ss ^ j dt j ((/3k + <y^^^)5un + <J^^}5ur^ dr, 

to dV(t) 

since the term containing dSur/dr can be integrated 
along the boundary and does not contribute. Vanishing 
total variation, 5Sb + SSg = 0, requires therefore 
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These are the known elastodynamic equations, together 
with the boundary conditions, now also considering sur- 
face tension. The velocity dependent term in the bound- 
ary conditions ensures momentum conservation. Notice 
that we had chosen for simplicity the particular case of 
a solid body inside the vacuum, and not for example in 
contact with a liquid or different solid phase. However, 
the calculations can easily be extended to these slightly 
more complicated case; then in Eqs. (|11() and l|12(l the 
difference [cr„i -|- Vnpui] between the values in the two 
phases appears instead [Tot- 



IV. THE CHEMICAL POTENTIAL 

In the previous section the dynamics of the mechanical 
variables, the elastic displacements, was studied, requir- 
ing that the variation of S with respect to Ui vanishes. 
The elastic body was already considered as moving or 
transforming, which leads to the appearance of the inter- 
face normal velocity in the elastic boundary conditions. 

The reason for this surface rearrangement is not yet 
specified and is not of mechanical origin. A specific ex- 
ample would be surface diffusion. The elastic motion as 
"mechanical" process describes the motion of individual 
atoms due to Newton's law; in contrast, the "chemical" 
processes are related to a (local) change of the number 
of atoms and consequently of the shape of the body. 

Demanding that the solid body is in equilibrium not 
only with respect to variations of the displacement vector 
but also with respect to this new degree of freedom, the 
interface position, leads to the requirement that SS = 
holds also for variations of the surface position in nor- 
mal direction by an interface shift Sn. Otherwise, small 
deviations from the condition 6S — initiate thermody- 
namical processes, driven by the chemical potential 
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From the bulk Lagrangian we gain 
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Notice that this expression contains in particular the dif- 
ference between potential and kinetic bulk energy. A 
naive guess would be that the chemical potential is the 
total energy density, containing thus the sum of both con- 
tributions. However, the above explanation shows that 
this is not true. A further discussion on this point is 
given in the next section. 

For the boundary contribution we assume here that 
the surface of consideration is free of external tractions. 
Then we get from Eq. HB4|I 
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where we used the rotation rule 
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and the expression for the total derivative 
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Therefore the total chemical potential, ii = iib + fJ-s, reads 
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Obviously, one immediately reproduces the known elasto- 
static case with surface energy, without inertial and sur- 
face tension effects ; notice that the chemical potential 
is then nothing else than the energy density at the sur- 
face. Inertial effects enter in an unexpected way, because 
instead of the sum of elastic and kinetic energy density 
their difference appears here. Surface tension gives rise 
to three contributions: At first, the term durr/dn re- 
flects the change of the strain Urr when the solid extends 
in normal direction, and Wtt-k expresses the elongation 
of a curved surface element. Without dynamical effects 
and external forces, the last term in the surface tension 
contribution vanishes, because, by virtue of the bound- 
ary condition (|12|l and Hooke's law, also Unr = 0. This 
full derivative appears only if both inertial contributions 
and surface tension are taken into account, leading to a 
coupling of both effects. 



V. THE DISSIPATION RATE 

An alternative approach to derive the expression for 
the chemical potential is to express the dissipation of 
energy due to morphological changes of the solid phase 
in terms of boundary integrals. Then the chemical po- 
tential gives the amount of energy change at the trac- 
tion free surfaces (per atom) if the shape of the solid is 



altered. It consists basically of three contributions: A 
change of the stressed volume reduces or increases the 
energy by the density il{Ub + T); also, a change of the 
shape modifies the boundary conditions and leads to an 
additional energy contribution — 2f2T. Finally, compli- 
cated "boundary" terms arising from surface tension and 
energy appear. 

To calculate the dissipation rate, we start with the bulk 
contribution, 



Eb^ j {T + Ub)dV. 
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Then we obtain 
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The boundary energy is 
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and its dissipation rate according to Eq. (|B3 
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Here the first integral is the contribution due to the 
change of the shape and constant fields, whereas the 
second one does the opposite, fixed shape and changing 
fields. We assume explicitly that the external forces are 
constant in time and that ?;„ = on the interface parts 
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where these forces act. From © follows for fixed shape 
(thus K = 0) in the second integral 

diir 
dr 

Using again Eqs. H14I) and (|15l) . we do not get a contri- 
bution from the full derivative, since it can be integrated 
along the whole (closed) surface. Thus we have 



dt 




Hence the total change of energy, E = Eb + Eg, is 
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The last integral gives the trivial contribution resulting 
from the external forces. We are mainly interested in 
local deformations of traction free surfaces {a^^^ — 0), 
leading to an alternative definition of the chemical po- 
tential. 
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with the same expression In particular, we clearly 

see that the density of the bulk Lagrangian instead of the 
total energy is relevant for the dynamics. 
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solid 
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FIG. 1: The Grinfeld instability. 



Here the boundary conditions are according to Eqs. (ini 
""nn — /3k and — 0. The latter one implies 
Unr = along the interface. This simplifies the general 
expression for the chemical potential, and we obtain 



fi = n[ ^a-^kUik -aK + 13 



dUrr 

dn 



The strategy is to solve the elastostatic problem for a 
slightly perturbed surface first, then to calculate the 
chemical potential, and finally to implement a dynam- 
ical process at this surface, driven by spatial variations 
of the chemical potential. 

Curvature in Cartesian coordinates is given by 

y" 



(1+^/2)3/2- 



For a two-dimensional plane-strain situation with isotro- 
pic elasticity we have 



<7xz ~ CTyz = 0, (Tzz = v{(^xx + yy) 



Strains are 



VI. GRINFELD INSTABILITY WITH SURFACE 
TENSION 

A straightforward application of the preceding results 
is the Grinfeld instability (ij ; similar to Q we also include 
surface tension here and neglect inertial effects. Usually, 
the dynamics of this instability, either driven by evapora- 
tion and condensation or by surface diffusion, is slow, and 
therefore kinetic terms need not be considered here. This 
is at least valid in the initial stage of the instability, but 
has to be corrected in the late regime when this omission 
leads to unphysical finite time singularities 0, Q . 

The geometrical situation is sketched in Fig. The 
horizontally applied uniaxial loading a^x = erg can be 
compressive or tensile. Our goal is to calculate the sta- 
bility of an originally flat surface. Therefore we study 
the evolution of a slightly perturbed interface, 

y{x, t) = A cos kx. 
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We use an Airy function U{x,y) to express stresses. 
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fulfilling the biharmonic compatibility condition A'^U = 
and the elastic boundary conditions anr = 0, (Jnn = /3k. 
We make the ansatz [111 



U = 



c^oy 



-iAiy+Bi)e''y cos kx+(A2y+B2)e^''y cos 2ka 



(20) 
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with Ai,Bi and ^2,-82 being of first and second order 
in the surface profile amplitude A respectively. Then we 
obtain from the boundary conditions up to second order 



Ai = -A(/3fc + ao), Bi=/?A, 



The total chemical potential at the interface is therefore 
to first order in A 



2E 



fcA cos fca 



-2ai 



E 



E 
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In the particular case of surface diffusion, 



(21) 



with a kinetic coefficient D, the amplitude evolution, A 
Aoexp(At), leads up to linear order to the spectrum 
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(22) 

assuming fc > 0. Obviously, the surface tension favors 
the Grinfeld instability. Apart from the Grinfeld length. 



La = 
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we have now an additional lengthscalc in the system, 

/3 



(23) 
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Consequently, the spectrum can be expressed as 
\L%ID = {kLG)'[l-ikLG) + l{kLG)^ 



Typically, Lp <^ Lq, since ctq ^ E, and thus the sur- 
face tension correction can be neglected. However, the 
self-stress due to surface tension in Eq. leads to an 
instability even without external loading do- 



VII. ENERGY CALCULATION 

In addition to the local chemical potential, also the 
total change of energy due to the morphological defor- 
mation of the surface can be calculated, leading to the 
same result. 



We compute the energy which is contained in a strip 
of the width 27r/fc. By symmetry, only even order con- 
tributions in A can appear. The total elastic energy, 

xo+2Tr/k y(x) 

Eb= j j Ubdydx 

x—xq y— — oo 

is infinite due to the zeroth order contribution (homoge- 
neous strain). We introduced a shift xq which is irrele- 
vant here since the energy content is the same in all strips 
of width of one period. The second order contribution in 
A is 



(2) _ _7rA_%_, _ 2^/_2 ,„2«2^ 
" ~ E 



The surface energy, 

Ec.^ j a\jl + y''^dx, 

X — Xq 

is to second order 
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Finally, the surface tension contribution 

Xf,-\-27r / k 

Ef3 = 



PUrrV 1 + y' dx 



(25) 



gives 



E^ = -^^^(1 - + 3ao). 

The expression H25(l can be presented in a more suggestive 
form by the means of Eq. , 

Uf3 = P[ur{xo + 2Tr/k,y{xo + 2TT/k)) - Urixo.yixo))] 

xq+2tt /k 

(SnUn dr. (26) 

X—Xq 

The first contribution vanishes for properly chosen offset 
Xq — 0. Also, it disappears for a closed boundary without 
corners. The second contribution is the work done by the 
total displacement against the surface tension (3k. 
Altogether, = E^^^ + E^J'^ + E^^ results in 



2ai 



A-v' 



E 



_fc^ii_i^(3ao + 2fc/3)). (27) 

The energy has the same structure as the spectrum (|22ll , 
therefore the prediction of the threshold of instability is 
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the same. We point out that the elastic constants appear 
only in the combination E/il — v^), since the scenario 
described here is plane strain. 

At a first glance, the description of the Grinfeld insta- 
bility using the chemical potential and the energy cal- 
culation seem to differ in a very important sense: The 
chemical potential is first order in the perturbation am- 
plitude A; in particular, it is sufhcient to calculate all 
elastic fields only up to first order. The use of energetic 
arguments always requires the calculation of second or- 
der terms. In principle, it needs the elastic fields also up 
to second order, because for example in the bulk energy, 

also terms like cr-°^M-^'' appear. This seems to contradict 
the former results using the chemical potential, which 
should contain the same physical information. In fact, it 
turns out that the second order terms in the elastic fields 
do not enter into the final energy terms; in other words, 
the same results can be obtained, setting A2 — B2 = Q va 
the Airy function above. The explanation is now rather 
simple for the surface tension integral: In the first rep- 
resentation (|25|l . it seems that Urr has to be evaluated 
up to second order. In the second representation (|25|l . 
assuming that the first term vanishes, the curvature is 
already of first order in A, and hence only the first or- 
der of Un is required. For the elastic bulk energy it is 
known that in the present case of linear elasticity there 
is no interaction between the stress caused by externally 
applied forces and internal locked-in stress by surface ten- 
sion |l2j | , and therefore only first order terms are needed 
here. 

Similar calculations of the total energy of a wavy sur- 
face were done in 0. However, the displacement fields 
without surface tension, as expressed by Eq. (5) there, are 
incorrect; the correct expressions can in principle easily 
be obtained from the Airy function l|20() . In particular, 
it would be necessary to include also third order har- 
monics cos Zkx to solve the elastic problem to the third 
order. Furthermore, the vertical displacement becomes 
correctly to first order (without surface tension) 



E~ 



o-qAcos kx -f 0(A^ 



in contrast to Eq. (5) in 2]. Additionally, in the expres- 
sion for the total energy of a system with surface tension, 
the work done by the surface tension against the (zeroth 
order) homogeneous displacement field was overlooked in 
0. For a closed contour, this energy contribution can be 
written as 



Up 



PKUn dr 



(28) 



dV 



with the total normal displacement u„ . For a non-closed 
surface also the corner contributions in Eq. H26I) have to 
be taken into account. 

This underlines the superiority of our formal approach 
to derive the total energy and the chemical potential, 
because intuitive guesses for the work done for example 
by surface tension are not required. 
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APPENDIX A: CONSERVATION OF 
MOMENTUM AND ANGULAR MOMENTUM 

Here we derive the basic conservation laws which result 
from continuous symmetries of the Lagrangian (^). We 
perform a variation of the fundamental field Ui, assuming 
prescribed evolution of the volume V{t) and dV{t). This 
means that a reshuffling of material by other processes 
than elastic displacement is possible; the total number of 
atoms does not need to be conserved. For simplicity, we 
assume the mass density p to be constant and no external 
forces a'l^ to be present. 

Using the expressions C3)-(01) and ||SJ), one gets 

6L ^ J {puiSui - aikSu^k) dV + j (inSun dr 

V[t) SV{t) 

piiiSui + ^^^5ui \ dV 
oxk J 



V[t) 



{PkSu„ - ainSui) dr 



9V{t) 



— J piiiSui dV, 
vit) 
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where the last step follows from the bulk equation HI()|I 
and boundary conditions (jlll II2|I . 

Obviously, the Lagrangian is invariant with respect 
to a rigid body translation, Sui = Si = const. ^ since 
Siii = and 6uik — 0. On the other hand, by Eq. IjAI |l . 
6L = EidPi/dt, with the total momentum 



pui dV. 



(A2) 



V(t) 



Therefore, the total momentum is conserved and the total 
surface tension force is zero. 

Similarly, we can consider rigid body rotations, 8ui = 
^ijkSjXk, with a constant axial vector e^; eijk is the to- 
tally antisymmetric tensor. These rotations do not affect 
strain, Suik = 0, and therefore again the Lagrangian 
is invariant, SL = 0, since also Siii — 0. By (|Aip . we 
have SL — Sjd£j/dt with the angular momentum 

J peijkXjUk dV. (A3) 

Vit) 

Thus the total angular momentum is conserved, too. 
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Notice that both conservation laws are not trivial, since 
the body can change its size and therefore its number of 
atoms: the normal velocity describes the local extension 
or shrinkage of the body, and is not specified here. 

An example is a solid droplet in its melt phase, for ex- 
ample in a temperature gradient. By crystallization the 
body can grow at one part of the interface and recede at 
another; this means that the droplet can wander, though 
all its "solid" atoms are at rest. In other words, the body 
moves, but its total momentum (|A2() vanishes, since the 
displacement is constant. 



APPENDIX B: AN INTEGRAL THEOREM 




We look at an energy functional of the form 



W[y{x,t),t] 



fdr 



V{t) 



FIG. 2: The interface at time t is depicted by the solid curve, 
the shifted interface ai t + 5t by the dashed curve. 



in Cartesian coordinates. Here the prime denotes the 
derivative with respect to x. We consider now the rate 
of change of the energy W: During a timestep 6t both 
the function value / and the interface position y change. 
Using the notations 



6y — TT- St, Sy' = Tr-TT St, 

^ dt ' ^ dxdt ' 



the rate of change is to first order in 5t 



5W = W[y{x,t + 5t),t + 5t]-W[y{x,t),t] 

I . ./ 

/2\l/2 



xo 

XI 



f V + (1 + v'^) 



(l+y'2)l/2 



5y 



oy oy ot 

(l + y'2)i/2_ 



oy ot 



>dx 



5y 



where we performed an integration by parts, assuming that the boundary contributions vanish; this is automatically 
fulfilled for a closed contour. The prime denotes the total spatial derivative (with fixed time i), i.e. 



dx 



dy ^ ^ dy' ^ 



Therefore 



5W 



dx dy dy' 



/2 // 

V y 



dy 



Sy+:;^{l + y'y/^5y 



dxdy' ' dydy'^ ' dy''^ 



y g _ r 



Curvature is given by 



and we have normal and tangential vectors 



1 



(l + iy'2)3/2^ 



(l + y'2)l/2 \^ 1 



1 



(l + y'2)i/2 yy' 
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Instead of Cartesian coordinates {x,y,y') (where y' is 
used to describe the slope), we use now local coordinates 
(n, r, 9). n is the normal, r the tangential coordinate and 
6 the inclination angle of the surface. With 6 — arctan?/', 
we get the following transformation rules: 



^ = (1 



d_ 

de 

d_ 

dn 
d_ 
dr 



y 



d_ 

dy'' 
y' d 



d 



[l + y'^yr^dx [l + y'^y/^dy' 



1 



d 



d 



[l + y'-^f/'^ dx (H-?/'2)i/2 
Comparing this with the above expression gives 



5W = 




df , 



drdd '*96'2 



6y 



^il + y'^)y^6t\dx. 



Orthogonality of normal and tangential variation requires 

6n^ = Sn^ + Suy = Sy 6ny . 

Furthermore, the normal vector is perpendicular to the 
interface {dx,dy), and therefore 



Suy = —Sn. 
" dr 



Combining the preceding expressions gives 



5n = 3^<5y, 
dr 



which, together with (|B2p . proves the statement ljBl|) 
above. 

Using this result, we finally obtain 



The next step is to show that a vertical shift 5y is related 
to a normal shift 5n according to 



Sy dx — Sndr — VnStdr, 



(Bl) 



where the last relation defines the interface normal ve- 
locity Vn- By substitution 



Sydx — Sy—dr 
dr 



(B2) 



with the arclength element dr^ — dx^ + dy^ (see Fig. \^ . 
Normal, tangential and vertical displacements are 



5t 



St,, 



Sn 



with Sn + St — Sy, thus 
St = 



Sn^ 
Sn., 



-Sllx 

Sy - Sny 



Sy 





Sy 



dW 

~dr 




V(t) 



dn ^ drdO de^\ ^ dt f 

(B3) 

The expression in brackets reflects the change of energy 
due to motion of the interface; in particular, it consists 
of a contribution caused by the motion of the interface in 
normal direction, the elongation of a curved length ele- 
ment and a tangential translational and rotational contri- 
bution; the last term is due to the change of the function 
/ in time. 

Similarly, we can derive the functional derivative 
SW/Sn, now at a fixed time t. It corresponds to w„ — 
S{xn), and thus 



sw df d'f d^f 

— ---/--T^--^- (B4) 



Sn dn 



drde 
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